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Abstract 

The internal s-electron polarization in the ground state of a Friedel- Anderson and 
a Kondo impurity is calculated. The Wilson basis with exponentially fine energies 
close to the Fermi level is used and expressed in terms of free electron waves with 
linear energy-momentum dispersion. It is found that in the singlet state the magnetic 
moment of the d-electron is screened by an s-electron cloud. The linear extension of 
the cloud is inversely proportional to the Kondo energy. When the singlet state is 
destroyed the polarization cloud is absent. 
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1 Introduction 



The properties of magnetic impurities in a metal is one of the most intensively studied 
problems in solid state physics. The work of Friedel pQ and Anderson [2] laid the foun- 
dation to understand why some transition-metal impurities form a local magnetic moment 
while others don't. Kondo [3] showed that multiple scattering of conduction electrons by a 
magnetic impurity yields a divergent contribution to the resistance in perturbation theory. 
Kondo's paper stimulated a large body of theoretical and experimental work which changed 
our understanding of d- and f-impurities completely (see for example [1], [5], [6], [7J, [8], [9], 
[TO] , [TTj . |12j . |13j). A large number of sophisticated methods were applied in the following 
three decades to better understand and solve the Kondo and Friedel-Anderson problems. 
In particular, it was shown that at zero temperature the Friedel-Anderson impurity is in 
a non- magnetic state. To name a few of these methods: scaling [TJ], renormalization |15j . 
[T6] , [TTj . [T8] Fermi-liquid theory [19], [20], slave-bosons (see for example [21J ) , large-spin 
limit [22], [23]. After decades of research exact solutions of the Kondo and Friedel-Anderson 
impurities were derived with help of the Bethe-ansatz [21], [25], [26], representing a magnif- 
icent theoretical achievement. The experimental and theoretical progress has been collected 
in a large number of review articles [7J, [8], [9], [10], [IT], [12], p3], [IS], [20], [21], [23], [21], 

m- m, m.- 

One of the most controversial aspects of the Kondo ground state is the so-called Kondo 
cloud within the radius where is called the Kondo length 
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(ksTx = Kondo energy, vf = Fermi velocity of the s-electrons) . 

The idea is to divide the ground state of a Kondo impurity into two parts with 
opposite d-spins. (By reversing all spins one can transform one component into the other 
one.) The proponents of the Kondo cloud argue that in each component there is an s-electron 
within the Kondo sphere which compensates the d-spin. This s-electron forms a singlet state 
with the d-spin. An important assumption of the Kondo-cloud proponents is that, above 
the Kondo temperature, the bond is broken and this screening cloud evaporates from the 
Kondo sphere. 

Already in the 1970's there were a number of theoretical papers with different predictions 
about the Kondo cloud which stimulated several experimental investigations. In most of the 
experiments nuclear magnetic resonance (NMR) was used, for example in Cu samples with 
dilute Fe-Kondo impurities. (For details see ref. [2B] which contains also an overview of the 
theoretical predictions at that time). By applying a magnetic field (small enough so that 
it does not destroy the Kondo singlet state) NMR was used to measure the electron spin 
polarization at shells of Cu atoms around the Fe impurity. The line shift (adjusted with 
the temperature dependent susceptibility of the Fe impurity) did not show any change when 
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the temperature crossed the Kondo temperature. This contradicted the concept that, in the 
Kondo ground state, the d-impurity was paired with and screened by an s-electron which 
evaporates above the Kondo temperature. 

In recent theoretical papers the argument is made that an NMR experiment cannot 
possibly detect the screening electron because of the large volume of the Kondo sphere. With 
a radius of the order of 0.1/xm (or larger) the Kondo sphere contains more than 2.0 x 10 8 
atoms. Then the Kondo cloud is only a faint Kondo fog. It is spread so thinly that the change 
of polarization felt by an individual Cu atom can not be detected in an NMR experiment. 

Actually the experimental side has not improved so far but a number of theoretical 
suggestions have been published since 2000 which propose to observe the Kondo cloud in 
sub- micron structures, in particular in connection with quantum dots [29], [30], [31], [32J, 
[33] . A number of these papers also contain numerical renormalization calculation to obtain 
the correlation functions between the impurity spin and the spin of the conduction electrons. 

In this paper I calculate quantitatively the screening cloud of a Friedel- Anderson and a 
Kondo impurity in the magnetic and the singlet state. 

The paper is organized as follows. In section 2 the theoretical background of this paper 
will be introduced, the FAIR method. In section 3 the numerical results for the spatial 
density of the s-electrons will be reported and discussed for several cases. First the spatial 
extent of the simple Friedel d-resonance with spinless electrons is discussed. Since the physics 
of the Friedel impurity is simple this consideration will give us confidence in the use of the 
Wilson wave functions. In the next step the electronic environment of the Friedel-Anderson 
(FA) impurity is investigated. The calculation is first performed for its magnetic state. The 
latter represents the magnetic building block of the singlet ground state of the FA-impurity 
which is investigated next. Last, but not least, the ground state of the Kondo Hamiltonian 
is analyzed. The subtle modification between the magnetic and the singlet ground state is 
discussed. Section 4 contains the conclusion. In the appendix details of the FAIR states and 
the Wilson states are presented, including the wave functions of the Wilson states and their 
density in real space. 



2 Theoretical Background 
2.1 The FAIR Method 

To investigate the question of a Kondo cloud I first consider a Friedel-Anderson impurity 
with spin 1/2 in a metal host. The underlying Hamiltonian consists of two Friedel resonance 
Hamiltonians (one for each spin) and a Coulomb term Und^n^i- 

{N N ~\ 

E 
+ E d S a d a + v sd{v)[dtc vc , + c+Xl \ + Un d,]n dA (2) 
u=l u=l ) 
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Here cJ, CT represent the creation operators of the Wilson s-electron state (which are discussed 
below), d\ represents the creation operator of the d-state at the impurity and n dtCT are the 
operators of the d-occupation. 

In the appendix I briefly describe the development of a very compact solution for the 
magnetic state and the singlet state. It uses two localized s-states a\ + and a\_ as artificial 
Friedel resonance states or FAIR states (Friedel Artificially Inserted Resonance states). The 
ground state solutions in terms of the FAIR basis have been studied in recent years for the 
Friedel Hamiltonian [31], [33], the magnetic and the singlet state of the Friedel- Anderson 
Hamiltonian [SB], [37] and the Kondo Hamiltonian [38J. The compact and explicit form 
makes these solutions ideal for calculation of spatial properties. The compact ground states 
which are used in this paper are given for the Friedel impurity by equ. (jSJ), the magnetic 
state by equ. d3J), the singlet state of the FA impurity by equ. (jlj) and the Kondo ground 
state by equ. 

The potentially magnetic state ^ms of the Friedel- Anderson Hamiltonian (which is the 
building block for the singlet state ^ss) has the following form 

n— 1 Ti—l 

^ms = [As^a^^l^ + A d:S dla* Q+ ^ + Ag^a^d* + A dtd d\d!\\ a*_ i a* +T $ (3) 

i=l i=l 

The magnetic state ^ ms has the same form as the mean field solution ^f m f- The only 
difference is that the state ^ms opens a wide playing field for optimization: (i) The FAIR 

states Oq + and aj_ can be individually optimized, each one defining a whole basis j a L j 

(which yields a Hamiltonian of the form in Fig. 10a in the appendix) and (ii) the coefficients 
A ai/3 can be optimized as well, fulfilling only the normalization condition A 2 S s + A\ s + A 2 S d + 
A\ d = 1. This yields a good treatment of the correlation effects. The optimization procedure 
is described in detail in ref. J3B]. The state is denoted as the (potentially) magnetic state 
$ms. It is expected to be a good solution above the Kondo temperature. 

The Coulomb interaction in the Friedel-Anderson impurity destroys the symmetry be- 
tween the occupation of spin up and down electrons and generates a magnetic moment (if U 
is large enough). However, this state with broken symmetry is degenerate. Its counterpart 
with reversed spins has the same energy, and at zero temperature the two form a new sym- 
metric state which has a different symmetry than original one. This is somewhat analogous 
to an atom in a harmonic potential. If some process transforms the harmonic potential into 
a double well potential the symmetry is broken. Although the atom has now an energy 
minimum on the left or right side its atomic wave function creates a new symmetric state 
by forming a symmetric superposition of left and right occupation. 

In a similar way I use the magnetic state as a building block for the singlet ground state. 
^ms together with its counterpart where all spins are reversed yield two states ^ms(TI) 
and ^ms (IT) • The singlet ground state is then given by 

#ss = ^ms(TI)±^ms(IT) 
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n— 1 n— 1 



^ S5 = [A SiS ag_ i ao +T + A diS d*ao +T + A^a^d*^ + A did d*d*] a*_^ a* +T $ (4) 



i=l i=l 
n—1 n—l 



+ 



v4 SiS ag +i ao_ T + A djS aQ H d* + A s , d d\a* Q _ } + A did d*d* 1 J a* +i a*_ T $ 





=i i=i 



(If all spin up states are moved to the left the plus sign applies). The bars on top of ^ms 
are supposed to point out that all the parameters such as al± and A a ^ have to be newly 
optimized. 

The resulting ground state is a very good approximation for the exact ground state. Its 
ground-state energy and the d-occupations for zero, one and two d-electrons at the impurity 
are of the same quality as the best numerical calculations in the field |39j . 

In the numerical construction of the FAIR ground state the states a} 0a is rotated in the 
iV-dimensional Hilbert space until the energy expectation value reaches its minimum. Each 
time one has to construct the other states of the basis j a L j so that they are orthonormal 

and non-interacting. This cannot be done with 10 23 electron states. Therefore we follow 
Wilson in constructing a reduced basis, in our case with generally iV = 50 Wilson states. 
The definition of the Wilson states, their properties and in particular their wave function in 
real space is discussed in the appendix. 



3 Numerical Results 
3.1 Friedel Impurity 

I start the analysis with the Friedel resonance of a spinless d-impurity. Since the Friedel 
Hamiltonian is a single particle Hamiltonian which can be solved exactly, one is here on 
familiar ground and the interpretation of the results will guide us in the less known territory 
of the Friedel- Anderson and the Kondo impurity. 

Given is the n-particle free electron state Wq = Yl c t^o i n terms of Wilson states. A 

(spinless) Friedel d-impurity is introduced at the position x = 0. This impurity has an s-d- 
interaction with the states ip v (x) = (c^ol^) • Due to the impurity the n-electron state 
is modified. As discussed in the appendix the exact ground state of the Friedel impurity can 
be written as 

V F = [Aa\ + Bd ] ) n«!$o (5) 

v ' i=l 

The composition of the FAIR state and the coefficients A, B are obtained by an iteration 
process in which the FAIR state is rotated in Hilbert space into its optimal orientation. 
Details of this iteration process and the construction of the full ground state have been 
described in previous papers [M] , [37] . 

In the next step we have to calculate the integrated electron density in the presence of 
the Friedel d-impurity at the origin. This calculation is discussed in the appendix. The 
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density of the state tp v (x) is given by 



p v (x) = \ip v (x)\ = 2 + '-dx 

TlX 

for v < N/2. In the evaluation we use N = 50 Wilson states. From the above equation one 
recognizes that the density p° u of the state ip v lies roughly in the interval between 2 U ~ 2 and 
2^+2 ^ n un j^ s of Air/2). Since for negative energies v takes the values from to (N/2 — 1) the 
different wave functions ip u (x) have very different spatial ranges and therefore very different 
densities, the lowest being of the order of 2~ 25 < 3 x 1CT 8 . Therefore one has to be careful 
in the summation of the different contributions. 

In the numerical evaluation we calculate (i) the s-electron density in the presence and 
absence of the d-impurity, (ii) form the net s-electron density p (x) as the difference, (hi) 
integrate the net s-electron density from x = to r where r is increased on an exponential 
scale, r = 2 l and plot the net integrated s-electron density versus I = log 2 (r) . This net 
integrated density q (r) = J Q r p(x)dx is plotted in Fig.l as a function of the distance from 
the impurity. The full triangles and circles give the integrated net s-electron density for the 

n—l n—1 

two components = Yl a l®o and ^/b = d) Yl aj^o °f the ground state, and the squares for 

i=0 i=l 

the full ground state typ- The curve for ^a increases from to 0.75 roughly between / = —2 
and +3. Around I — 23 it drops back to zero. It has to return to zero because the total 

n-l 

number of occupied states in ^a is the same as before in the state = EI c l®o, i-e. equal 

i=0 

to N/2. The equivalent distance r = 2 23 corresponds roughly to the maximum range of the 
wave functions. In a way it can be considered as the border or surface of the theoretical 
sample. Therefore the interpretation of this curve in Fig.l is that in the state a fraction 
of 0.75 electrons has moved from the surface towards the d-impurity. There it surrounds 
the empty d-state in the spatial region between r = 2~ 2 and 2 3 . (The component ^a has 
no d-electron). The behavior of the integrated net density beyond about I w 20 is a surface 
effect and of no interest for the present investigation. If one increases the number of cells 
from iV to N' then this part of the curve moves by (N' — N) /2 to the right while the left 
part of the curve remains unchanged. 

The full circles in Fig.l give the change of the net integrated density of s-electrons for 
the state ^ B as a function of distance. At large distances (r w 2 23 ) it approaches the value 
— 1 because this state has only (N/2 — 1) s-electrons. Between the distances / = —2 and +3 
the curve for tyg assumes the value —.25. If one adds to the integrated s-electron density 
the one d-electron at the origin then one obtains the dotted curve which starts on the left 
at one and approaches zero at large distances. It is interesting to notice that in the range 
5 < I < 20 the total integrated density of d- and s-electrons is equal to 0.75 for both wave 
functions ^a and ^b- This means both components have accumulated the same charge close 
to the impurity. As discussed below this corresponds to a Friedel phase shift of 5 = 0.757T. 

In this specific example the parameters of | 2 = 0.1 and Ed = —.135 were chosen so 
that the occupation numbers take the values A 2 = .25 and B 2 = .75. Then the net integrated 
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density of ^a times the occupation number A 2 yields 0.75 * 0.25 « 0.188 and just cancels 
the corresponding product of 0.75*(— .25) for the state ^b- As a result the total integrated 
density of s-electrons as a function of distance from the impurity vanishes as the curve for 
typ shows (full squares). There is no change of s-electron density about the impurity. On 
the other hand the curve approaches the value of —.75 for large distance. This is the charge 
which has been moved from the surface into the d-state. This charge is predicted by the 
Friedel sum rule. The occupation B 2 = 0.75 of the state ^b means that the d-occupation is 
0.75, i.e. the phase shift is S — .75n and therefore the Friedel charge is 5/ir = 0.75. 
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Fig.l: The net integrated density f^p (x) dx is plotted versus the logarithm 
of the distance / = log 2 (r) from the impurity for the components ^a, 
and the full wave function ^/ F . For the dotted line the additional one 
d-electron has been added. 

The Friedel resonance is also ideal testing ground for the spatial extension of the wave 
function. For that purpose the d-state is moved to the Fermi level (E^ = 0). The half 
width of the resonance is given by A = n \ V® d f po (po=density of states). This resonance 
width introduces an additional length scale into the problem. In the resonance the d-state 
hybridizes with the s-electrons in this energy range of 2A.This corresponds to a wave number 
range of k& ~ 27r |V^| 2 po an d therefore a smearing of the wave function over a range xa oc 

1/ (^2n |V^| 2 po) (Vgd 1S measured in unit of Ef and xa in units A^/2). 

In the numerical evaluation I change the value of jV^ | 2 in steps of two from 0.005 to 0.16. 
Indeed one recognizes in Fig. 2 that with decreasing \V® d f the increase of the integrated net 
density is stretched over an increasing length scale. Since the horizontal separation of the 
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points by 1 corresponds to a factor of two in the length scale each curve is broader by a 
factor two than the previous one. 



Friedel state 

0.6- 




ra-l 

Fig. 2: The net integrated density J Q r p (x) dx for the state ^a= Yl a l^o for different values of 

i=0 

\Vg d \ while the d- resonance lies at the Fermi level, E d = 0. The width of the wave 
function increases proportional to 1/ir \V® d \ . 

3.2 The Magnetic State 

The magnetic state ^ms is the building block of the singlet state. Its multi-electron state 
is built from four Slater states and shown in equ. I calculate the net integrated density 
of spin up and down, as well as total density and spin polarization in the vicinity of the 
impurity when it is in the magnetic state. For the present purpose I choose the parameters 
E d = 0.5, \V s ° d \ = 0.04, U = 1 and iV = 50. This yields a well developed magnetic moment 
of fi = 0.93/iB- The occupation of the different components is A? s s = 0.0294, A? sd = 0.0057, 
A\ s = 0.9355 and A\ d = 0.0294. 

In Fig. 3 the net densities are plotted. As discussed above, the range beyond 2 20 corre- 
sponds to the rim or surface of the sample and is of no interest for the density distribution 
around the impurity. One recognizes that there is only a negligible polarization of the elec- 
tron gas in the vicinity of the impurity. The important result of Fig. 3 is that there is no 
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polarization cloud around the magnetic state of the impurity. 
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Fig.3: The net integrated density p(x) dx of the s-electron within a distance r 
from the impurity for spin up and down, as well as total density and spin 
polarization. The magnetic moment of the impurity is 0.93^. 



3.3 The Singlet State of the Friedel- Anderson Impurity 

The author derived the compact singlet ground state in ref. [37]. Its form is given in equ. (TjJ. 
The singlet state does not have a net spin polarization because of the symmetry between spin 
up and down. However, the individual components may possess a spin polarization cloud. 
To investigate this question one has first to decide how one subdivides the symmetric state 
of the impurity. Since there are four Slater states which have either zero or two d-electrons 
one cannot divide the full state into two parts in which either the d^- or the <ij-state are 

occupied. A reasonable division would be into the states ^ ms (Tl) an d ^ ms (IT)- However, 
one may object that both states contain a component in which, for example, the (i-f-state is 
singly occupied. Nevertheless I use the components ^ms(M) as the magnetic component 
because in the singlet state the coefficients A SjS , ..A dd and A s>s , ..A djd have the same sign 
while in the triplet state they have opposite signs. (Their absolute values are in both cases 
pair- wise identical). 

In the following plots in Fig.4 I choose the state ^ms (Tl) with the majority d-spin 
pointing up. The same parameters as for the magnetic state are used: E d = —0.5, \V® d \ 2 = 
0.04, U = 1. This yields the following coefficients: A 2 SS = 0.0146, A 2 sd = 0.0028, A 2 ds = 0.4629 
and A 2 dd = 0.0146. They are about half the values of the magnetic state. Therefore on a 
superficial glance it appears that ^ ss w (l/v 7 ^) [^ms (Tl) + *ms (IT)] • 
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In Fig. 4 the integrated densities of spin up and down electrons, their sum and difference 
(the polarization) are plotted versus the distance from the magnetic impurity (on a logarith- 
mic scale). One recognizes that now one has considerable contributions to the integrated 
net densities of both spins. The polarization of the two contributions is no longer zero but 
reaches a value of —0.46 at a distance of r « 2 116 . Since the magnetic state with net d-spin 
up has only a weight of about 1/2 it contributes an effective dj-moment of 0.93/2 « 0.46. 
Therefore this d-spin is well compensated by the polarization of the s-electron background. 
The difference with the pure magnetic state is particularly striking. We observe a screening 
polarization cloud of s-electrons about the impurity within the range of r « 2 116 = 3.1 x 10 3 . 
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Fig. 4: The net integrated density J Q r p (x) dx within a distance r = 2 l 
from the d-spin up component of the impurity. Shown are the spin up, spin 
down components as well as the total density and the polarization. The d^-spin 
of 0.93/2 is screened by 0.46 s-electrons within the range of r w 2 1L6 . 

The range of the screening cloud depends on the strength of \V® d | 2 . In Fig. 5 a similar 
plot as in Fig.4 is performed for an s-d-coupling of |V^| = 0.1. This time the polarization 
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cloud extends only over a distance of r » 2° » 84. 
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Fig. 5: The net integrated density J Q r p (x) dx within a distance r = 2 l 
form the d-spin up component of the impurity. Shown are the spin up, spin 
down components as well as the total density and the polarization. The d^-spin 
polarization cloud extends for \V® d | 2 = 0.1 over a range of r w 2 6,4 . 



3.4 The FAIR states 

It is remarkable how different the polarizations of the s-electrons in the magnetic state ^ ms 
and the magnetic component of the singlet state ^ ms are. This is particularly surprising 
since the amplitudes A a $ = A a ^ in the singlet state are roughly equal to A aj p / \/2. However, 

in the singlet state one has a finite coupling between ^ms and ^ms- This affects the 
composition of the FAIR states a) Q+ and al_. The two states and aj_ contain the 
whole information about the many electron states ^ms or ^ss- When known the 

whole bases ja| + j and |<4_ j and the coefficients A n ,a can be reconstructed. Therefore it is 



worthwhile to compare the two FAIR states for the magnetic and the singlet state. 

For the magnetic state the coefficients «q ± of the states aj ± = J2 u cto±ct are plotted in 
Fig. 6 versus the cell number u, which is a measure of the energy. While a^ Q+ is essentially 
concentrated at positive energies the coefficients of Oq_ have their main weight at negative 
energies. The two are mirror images with respect to the energy zero. 
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Fig.6: The coefficients of the FAIR states a Q+ and a _ are plotted in terms 
of the states cj, for the magnetic state The full and the dashed curves 

represent aj + and a\_ respectively. 

i 1 2 

Fig. 6b: The density distribution of the p v = poi /^^ as a function of as a 
function of v. 

In Fig. 6b the quantity Pj,=|o:q| 2 /A u is plotted as a function v, where A u = (C„+i — (u) is 
the width of the energy cell <£ u of the Wilson state ip u (x) . The quantity p u is a function of 
the energy p v = p (C) , and p (() d( represents the weight of the original basis states <fk (x) 
in the energy window (£, £ + d() to the FAIR state aj. This weight is plotted in Fig. 6b for 
the magnetic state as a function of v. It would be more natural to plot p (() as a function of 
the energy. But for v close to N/2 the width of the energy cells <£ u is less than 1CT 6 and any 
dependence of p (() on the energy cannot be resolved on a linear scale. The probability p v 
increases close to the Fermi energy. Again the contributions of aj + and a\_ resemble mirror 
images. 

In Fig. 7a and Fig. 7b the corresponding plots are shown for the magnetic component of 
the singlet state. While the overall shape of the coefficients in Fig. 7a resembles that in Fig.6 

I 1 2 

the weight p u ± = |«o±| I close to the Fermi surface is very different for the singlet state 
and the magnetic state. While the two probabilities are mirror images for the magnetic 
case and reach a maximum value of about 10 one observes in the singlet state a maximum 
of about 400 and the weights p u ± in a\ + and aj„ are essentially identical and not mirror 
images. The magnetic component of the singlet state is in a subtle way different from the 
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magnetic state. 




Fig. 7a: The coefficients a^ ± of the FAIR states a 0+ and a Q _ in terms of 
the Wilson states cf, of the magnetic component ^ms of the singlet state. 

I 1 2 

Fig. 7b: The density distribution of the p u = «o± /A y as a function of 
the cell number v. 

It is interesting to note that the spatial densities of the two FAIR states a\ + and aj_ 
(after averaging over Friedel oscillations) are identical with an accuracy of four digits. This 
applies for the magnetic and the singlet state. However, their energy expectation values 
are essentially opposite equal. The values for ^<2q + $o \H\ c^+^oj are rather similar in the 
magnetic state (~ +0.288) and the singlet ground state (+0.2865) of the FA impurity. 



3.5 Kondo Impurity 



If one increases the Coulomb potential in the Friedel-Anderson impurity (and lowers the d- 
state energy Ed) then the coefficients A SyS and A^a converge towards zero. This is the limit 
of the Kondo impurity. By means of the Schrieffer-Woulff transformation one can express 
the s-d-interaction as an exchange interaction of the form 2 Js • S. In this limit the compact 
ground state takes the form 



K 



A s4 a\ + ^d\ + A djS d\a\ ) _ i \[ aj +T f] a ]-T $0 



n-l 



n-l 



(6) 



+ 



i=i 

n-l 



3=1 
n-l 



J i=l 



Since equ. (jSJ) is a reduced form of equ. (pEJ) the calculation of the integrated net density is 
analogous to the previous calculation. The Kondo singlet state is composed of two magnetic 
components represented by the first and the second terms. Note that the effective moment 
of these magnetic building blocks is less than \\xb- 

We investigate the Kondo ground state for the exchange parameter J = 0.1. In Fig. 8 the 
net integrated density is plotted versus the logarithm of the distance r = 2 l for the magnetic 
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component with the d-spin up. One recognizes that a spin polarized "cloud" extends up to 
a distance of r = 2 14 7 from the Kondo impurity. 
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Fig. 8: The net integrated density J Q r p (x) dx within a distance r = 2 l 

from the d-spin up component of the impurity. Shown are the spin up, spin 

down components as well as the total density and the polarization. 

Next we want to study the extent of the Kondo cloud as a function of the Kondo energy. 
In a recent investigation we obtained for the coupling of J = 0.1 a relaxed singlet-triplet 
excitation energy of AE st (0.1) = 2.35 x 10 -5 . By reducing the coupling to J = 0.08 we 
observed an AE st (0.08) = 1.37 x 10~ 6 . The ratio of the two energies is equal to 17.1. 

Fig. 9 shows the integrated net s-electron density for a Kondo impurity with J = 0.01. 
The two polarization curves (stars) in Fig. 8 and Fig. 9 are very nicely parallel. Therefore the 
relative shift of the two curves can be determined very accurately. It is equal to 4.1 in terms 
of the log 2 -scale, which corresponds to a length ratio of 2 4 ' 1 = 17. 1. This is a very good 
agreement with the ratio of 17.1 of the two Kondo energies. Therefore we can confirm that 
the extent of the screening cloud scales with the Kondo energy. Kondo energy and extension 
in real space are reciprocal. 
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Fig.9: The net integrated density f Q p (x) dx within a distance r = 2 l 
from the magnetic component of the impurity. Shown are the spin up, spin 
down components as well as the total density and the polarization. The 
polarization extends over a distance of r = 2 18 8 s; 4.6 x 10 5 [Xf/2] . 

If in Fig.9 we extrapolate the range of the Kondo cloud as the intersection between the 
linear decrease and the saturation at the right then we obtain an intersection of / = 18.8 
or a range of 2 18,8 (in units of Ai?/2). In the reduced length units of Ai?/2 and the energy 
units of Ep the Kondo length equ. ([1]) has the values £x = 2/ (ttek) where Ek is the Kondo 
energy in units of Ep. For the parameters of Fig.9 (J = 0.08) the Kondo energy has the 
value e k = AE st = 1.37 x 10~ 6 . This yields £ K = 4. 6 x 10 5 « 2 18 8 . The perfect agreement 
with the numerical extrapolation is probably lucky because the equation ([T]) for the Kondo 
length is intended as an order of magnitude. 



4 Conclusion 

A recently developed compact solution for the singlet state of the Friedel- Anderson and the 
Kondo impurity is used to investigate the old question of a Kondo cloud in the Kondo ground 
state. Wilson's states with an exponentially decreasing frame of energy cells towards the 
Fermi level and a linear dispersion relation are expressed in free electron waves integrated 
over the width of their energy cells. A rotational iteration process (which has been described 
in previous papers) is used to obtain the optimal ground state. The singlet state of the 
Friedel- Anderson impurity consists of eight Slater states. They are divided into two groups 
of four, where each group represents one magnetic component of the singlet state. Two 
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criteria yield the same splitting of the eight Slater states: (i) Each group transforms into the 
other one by inverting all s- and d-spins, (ii) The two groups combine symmetrically in the 
singlet and anti-symmetrically in the triplet state. 

After identifying the magnetic components, the real space wave function is evaluated. 
Since the spatial range of the Wilson states varies in the range between a Fermi wave length 
Xf and 2 N ^ 2 \p (N « 50 is the number of Wilson states used in the evaluation) one has to deal 
with electron densities which vary be a factor of 2 25 « 3 x 10 7 . Therefore it is favorable to 
integrate over the net s-electron density from the impurity outwards. The resulting curves 
yield the spatial range over which the spin densities of the s-electrons contribute to the 
polarization about the impurity. 

If one wants to obtain the polarization about the impurity then one has to calculate 
the difference of spin up and down s-electron densities with and without the impurity. The 
simple and exact solution of the Friedel resonance for spinless electrons is used to test the 
method and interpret the results. It turns out that for each Slater state the behavior of 
the (integrated) density and polarization can be divided into two regions. For r > 2 N / 2 the 
system follows the sum rule which is imposed by the number of s- and d-states in the multi- 
electron Slater state. For r < 2 N I 2 one obtains the physical densities and polarizations. The 
behavior in this range does not change when N is increased. If the d-state in the Friedel 
resonance is set at the Fermi level then one observes in real space that the range of the 
s-electron cloud is proportional to 1/ |V^| 2 . 

For the magnetic state m s of the Friedel- Anderson impurity one observes essentially no 
magnetic s-polarization in the vicinity of the d-impurity. 

For the magnetic component ^ms of the singlet state in the Friedel- Anderson impu- 
rity one observes an s-polarization cloud which screens the spin (magnetic moment) of the 
d-electron. The range of this polarization cloud is investigated in detail for the Kondo im- 
purity. The range is inversely proportional to the Kondo energy. The latter was obtained 
in a previous investigation as the energy difference between the singlet state and the re- 
laxed triplet state energy. The absolute value of the range of the polarization cloud agrees 
surprisingly well with the Kondo length. 

The different screening behavior in the magnetic state ^ ms and the singlet state ^ ss is 
due to subtle differences in the composition of the FAIR states aj + and a^_. In the singlet 
state the FAIR states have a much larger weight (of the original basis states ipk) very close 
to the Fermi energy. It is remarkable that the two FAIR states aj + and Oq_, which are far 
apart in their energy, have essentially an identical density in real space (after averaging over 
the Friedel oscillations). 

In the present paper a very simple energy band and dispersion relation is used in analogy 
to Wilson's work. This simplifies the numerical evaluation dramatically. However, the spatial 
dependence can be evaluated for an arbitrary s-band with an energy-dependent density of 
states and s-d-interaction. The author generalized the definition of the Wilson states [10] 
for this case. The numerical work would, however, be much more extensive. (The spin of the 
impurity is still restricted to 1/2). 

It should be emphasized that the extent of the Kondo cloud calculated here applies only 
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when the mean free path of the conduction electrons is larger than the Kondo length. Also 
the surface of the host should be at least this distance away from the impurity. This does, 
however, not mean that for smaller sample size the Kondo effect is suppressed. As long as 
there are sufficiently many electron states within the Kondo resonance at the Fermi level 
which couple to the impurity one obtains the full Kondo effect [31] ■ 

Finally it should be noted that the extent of the electron density in real space is a 
detector for a resonance in energy. The spatial extension £ and the resonance width A are 
reciprocal and given by the relation £A « Hvp. The Wilson states act here as a magnifying 
glass close to the Fermi energy because their energy width A u decreases as 2~ v towards 
the Fermi level. This is particularly demonstrated for the Friedel impurity in Fig. 2. The 
extent of the net integrated density is only observed when the energy of the d-state lies at 
the Fermi level. Otherwise it cannot be detected. For example with Ed = —0.5 the the 
resolution of the Wilson states is only between 0.5 and 0.25 which is not sufficient to detect 
the Friedel resonance. One has to split the Wilson states close the resonance to make it 
visible. Therefore the Wilson states can be used as a quasi-experimental spectroscope. 

A The FAIR Method 

The Hamiltonian of the Friedel- Anderson impurity is given in equ. (J2J). One obtains 
the mean-field Hamiltonian from equ.(T5]) by replacing rid+rid- —>n d +(nd-) +(nd+)nd- 
— (nd+) {rid-). After adjusting (rid+) and (rid-) self-consistently one obtains two Friedel 
resonance Hamiltonians with a spin-dependent energy of the c? cr -state: E dja = E d + U (n d - a ). 

{N N 1 

u=l u=l ) 

The mean-field wave function is a product of two Friedel ground states for spin up and down 

It has been shown [33] , [55] that the exact n-particle ground state of the Friedel Hamil- 
tonian can be expressed as by the sum of two Slater states, in which either the cf'-state or 

t n_1 t 

a state a is multiplied with the same (n — l)-s-electron state Yl a }®o (see equ. (J5])). The 

i=i 

state aj = J2u a o^ c l ^ s a localized state whose components a 0)U can be obtained analytically 
[35] or by variation [33]. 
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Fig. 10a: The matrix of the Friedel Hamiltonian in the basis j a ! j > $ ■ 
Fig. 10b: The two basis vectors a , d) are rotated into x^ = \Aa\y + Bd) 



and y^ = (Bclq — Ad)j . For the optimal aj the matrix is block-diagonal. 

In Fig. 10a the matrix of the Friedel Hamiltonian in the basis | a l| > ^ is shown. Only the 

diagonal elements and the elements along the edges are non-zero. The inner (N — 1) x 
(N — l)-sub-matrix is diagonalized. This way the state determines uniquely the full basis 

| a l|- O ne recognizes that the state represents an artificially inserted Friedel resonance 

state, with similar properties as the state <$ . Therefore I call a\ a "Friedel artificially inserted 
resonance" state or FAIR-state. The use of the FAIR-states is at the heart of my approach 
to the FA- and Kondo impurity problem. Therefore I call this approach the FAIR method. 

The structure of the Hamiltonian in Fig. 10a can be obtained for an arbitrary state a\. 
There is, however, one optimal state which yields the exact n-electron eigenstate p 
and for which the energy expectation value (typ \Hp \ typ) has a minimum. With this state 
the two basis vectors aj, d) can be rotated within their plane into the orthogonal states 

= \ Aclq + Bd^j and = ^-BaJ — Adyj . For the optimized aj the resulting Hamiltonian 

is block-diagonal. This is shown in Fig. 10b. The upper left block contains the occupied 
states. Since there is no matrix-element between the occupied and the empty block the 

product state x^ala^-^n-i^o = (^ a o + Bd)\ Yl a l^o is an eigenstate. 

For the two Friedel states in the mean-field wave function I use the form of equ. (jSJ) and 
obtain for the mean-field solution 



rnf 



B + d\ 



n-1 

n 

i=l 



A A-i + B-d\ 



n-1 

n 

i=l 



a 



(7) 
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where j and | a I-| are two (different) bases of the ^-dimensional Hilbert space. This 
solution can be rewritten as equation (jSJ). 

In the mean-field solution ^ m f the coefficients A a ^ are restricted by two conditions 
A\ + B\ = 1 (here A SyS = A + A_ etc). Therefore this state does not describe well the 
correlation effects. 

In contrast the state ([3]) opens a wide playing field for improving the solution: (i) The 
FAIR states aj + and a\_ can be individually optimized, each one defining a whole basis 

f a l±} (which yields a Hamiltonian of the form in Fig. 10a) and (ii) the coefficients A a ^ can 
be optimized as well fulfilling only the normalization condition A 2 S s + A\ s + A 2 S d + A\ d = 1. 
This yields a much better treatment of the correlation effects. The resulting state is denoted 
as the (potentially) magnetic state ^ms- The magnetic state ^ms has the same form as 
the mean field solution \I/ m /; the only difference is that its components are optimized for 
the Friedel-Anderson Hamiltonian. The optimization procedure is described in detail in ref. 

The magnetic state ^ms is used as the building block for the singlet state. It form is 
given in equ. (j4]). 



B Wilson's states 

Wilson considered an s-band with constant density of states and the Fermi energy in the 
center of the band. By measuring the energy from the Fermi level and dividing all energies by 
the Fermi energy Wilson obtained a band ranging from —1 to +1. To treat the electrons close 
to the Fermi level at £ = as accurately as possible he divided the energy interval (—1 : 0) 
at energies of —1/2, —1/4, —1/8, .. i.e. ( u = —\j1 v . This yield energy cells t v ranging from 
-1/2" to -l/2^ +1 with the width A„ = ( u+1 - C* = l/2" +1 . 

Wilson rearranged the quasi-continuous original electron states <fk (x) in such a way that 
only one state within each cell <L V had a finite interaction with the impurity. Assuming that 
the interaction of the original electron states ipk (x) with the impurity is /^-independent this 
interacting state in £ v had the form 

(x) = Ee,^fc (x)/y r Z^, 

where Z u is the total number of states ipt (x) in the cell <£ u (Z v = Z (Cv+i — Cv) /2, Z is the 
total number of states in the band). There are {Z v — 1) additional linear combinations of the 
states (fk in the cell but they have zero interaction with the impurity and were ignored 
by Wilson as they will be within this paper. 

The interaction strength of the original basis states ipk (x) with the d-impurity is assumed 
to be a constant, v S( i- Then the interaction between the d-state and the Wilson states i[) u (x) 
is given by V sd (u) = KdV (C+i - &0 / 2 where \ V sdt = E fc \ v sd\ 2 = \ V sd (^)| 2 - 
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B.l The wave function of Wilson's states in real space 

For the discussion of the wave functions in real space one has to look a bit closer. We assume 
a linear dispersion relation between energy and momentum, a constant density of states and 
a constant amplitude at x — 0. Then the interaction between the original basis states <fk (x) 
and the d-impurity will be constant for all k-states. These assumption are the same as in 
Wilson's treatment of the Kondo impurity We define the wave functions (x) is such a 
way that the results apply for the impurity problem in one, two and three dimensions. 

B.l.l One-dimensional case 

Let us start with the one dimensional problem. Here we have the impurity at the position zero 
and the conduction electrons are located in the range between and L. The wave functions 
(fk (x) have the form ip k (x) = a/2/Lcos (kx). There is another set of eigenstates ip k (x) = 
A/2/Lsin (kx). These states don't interact with the impurity at the origin. Therefore they 
don't have any bearing on the impurity problem. 

B.1.2 Three-dimensional case 

In three dimensions the free electron states can be expressed as ip k (x) = Y™ (8, <p) ji (kx) 
where Y™ is a spherical harmonics and I, m are the angular momentum and magnetic quan- 
tum numbers, ji (kx) is a spherical Bessel function. Its long range behavior is given by 

(l/r\ / 2nLj sin (kx — In/ 2). Only the states with the same I as the impurity couple to the 

impurity. All the other states for different I belong to the group of inert states ip k (x) . 

If one calculates the density of the wave function, integrates in the three-dimensional 
case over the spherical surface 4nr 2 and averages over short range (Friedel) oscillations then 
one obtains in the 



• one-dimensional case: (2/L) cos 2 (kx) = 2/L 

• three-dimensional case: p k (x) = (2/L) sin 2 (kx — Iti/2) = 2/L 

In both cases one obtains essentially the same density. Therefore it is sufficient to use 
the one- dimensional approach for calculating the density of a Kondo cloud. It is equivalent 
to the 3-dimensional case integrated over the spherical surface. 

B.1.3 The wave functions in one dimension 

While the energy is measured in units of the Fermi energy the momentum will be measured 
in units of the Fermi wave number. We assume a linear dispersion relation for < k < 2 
with 

C = (*-i) 
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The (almost) continuous states ipk are given as 



(x) = y^ cos (nkx) 



where L is the length of the one-dimensional box. Since k is measured in units of hp the 
coordinate x gives the position in units of Xf/2 where A^ is the Fermi wave number. The 
boundary condition cos (irkL) = yields k — (A + 1/2) / L. (The maximal value of A is 2L, 
since k is dimensionless then L is also dimensionless) . Therefore we have Z = 2L states in 
the full band of width 2. 

To obtain the Wilson state we have to sum the states ipk (nkx) over all states within an 
energy cell. If the cell ranges from (£„ : C^+i) corresponding to a /c-range (1 + ( u ) < k < 
(1 + (u+i) then we represent all the states in this energy interval by 

From Z u = L (C^+i — Cu) states we have (according to Wilson) constructed one state ip v (x) 
which couples to the impurity. The same procedure yields [Z v — 1) additional linear com- 
binations of ifik (x) which do not couple with the impurity at the origin. We denote these 

states as <fk (x)- They are as inert to the impurity as the states <fk (x) and will be included 
in the quasi- vacuum. 

The wave function of the state of the states cj, has the form for v < N/2 

, , , 2y^ sin (g^g) f 7Tx{2 + Cu + Ul] 

Wv (X) = = — COS — 

- C,) ™ V 2 

For the exponential energy scale this yields for v < N/2 — 1 

1> v (x) = 2V¥^ Sin cos (vx (l - ^) ) (8) 

Similarly one obtains for in the positive energy range 

1> N + V (x) = 2V2=^ Sin cos (vrx (l + JL) ) 

The two wave functions ipN/2-1 and ipN/2 are special because their /c-range is the same 
as their neighbors ipN/2-2 an d N/2+1 All four states close to the Fermi level have the same 
/c-range of 2~ Ar / 2 ~ 1 . One has to pay special attention to this complication. 
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B.2 Density of the Wilson states in real space 

The density of the wave function ip v (x) is given by 

,2 ( 7TX-(C^ + l-Cy) 



sin' 

2 



hMz)l =77 : — T2 cos 



TlX (2 + Cu + Cu+l) 



The density of a single state ip u (x) is given by the square of the function if)„ (x) in equ. 
(JHJ). This density has a fast oscillating contribution which yields the Friedel oscillations. We 
are here interested in the density on a much larger scale and average over the fast oscillation 
(which has a period of 1 in units of \p/2). Then we obtain for the Wilson states 



, sin* 



(jCX 2 i/+2 ) 



|VV (x)\ 2 = 2"+ 3 - V '~ 2V+ ~' J dx (9) 



TCX 



In the numerical calculation we will use (most of the time) N = 50 Wilson states. From 
equ. ([9]) one recognizes that the density p° v of a state ij) v lies roughly in the interval between 
2 U ~ 2 and 2 U+2 (in units of Xp/2). Since for negative energies v takes the values from to 
(N/2 — 1) the different wave functions ip u (x) have very different spatial ranges and therefore 
very different densities, the lowest being of the order of 2 -25 < 3 x 10 -8 . This means that 
it is not useful to calculate the density as a function of x because this density varies over a 
range of 2 25 . Instead, we use the integrated density, integrated from to r. 

gO(r)= [ \^(x)\ 2 dx = 2^ ^ Sln2 dx 



Jo Jo nx 

r rf 

f 2^+ 2 sin (iru) , 

= 2 / -rV d « 

Jo (iru) 

One realizes that a single integral yields the integrated density for (almost) all wave function 
ipv (x). The state i/jq (ttx) lies roughly in the range between 2~ 2 < r < 2 2 , i.e. the integrated 
density q u (r) = j Q r \tp (ttx)\ 2 dx increases in this range from 10% to 90%. Therefore the 
states ipv (x) and ipN-i-v (itx) lie in the range between 2 U ~ 2 < r < 2 U+2 . For a total of 
N = 50 Wilson states the maximum range of the wave functions is roughly 2 N I 2 = 2 25 . 
We may define as a ruler a linear array I (s) where s is an integer, (—N/2 < s < N) as 

i( S )=2 r^p^-du 

Jo (iru) 

Then the integrated density of the state i/j u in the range from to 2 l is given by 

<7° (2') = f |VV (x)\ 2 dx = 2 ^p^fidu = I (I - v - 2) 

Jo Jo (iru) 

Then g° (2') gives the integrated density on an exponential scale. 
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B.2.1 Interference terms in the density 



The Wilson states ip u (x) or cj, represent the free electron states in the impurity problem. 
After the interaction we express the ground state in terms of new states a\ = J2pZo a i c i- 
Its integrated density is given by 



Pi (2 1 ) = f YH=o<^ (x) 
Jo 



2 

dx 



The quadratic terms can be evaluated with the same ruler / (s) as before. But this 
time one has in addition interference terms ip u (x) (x) . These terms depend on two 
parameters, v and A. So one needs for each A a different ruler. Furthermore the interference 
terms depend on the sub-bands of ip u (x) and Vv+a (%)■ If both lie either in the negative 
energy sub-band (v, v + A < N/2) or in the positive sub-band (v, v — A > N/2) then one 
obtains one set of rulers I (z/, A) and if they lie in opposite sub-bands then one obtains 
another set of rulers Ji (z/, A). As an example one obtains 



J (S, \)=2V¥ V { 2X> COS 37TO 1 - - 

J0 (7TU) V V 1 



du 



For 1\ (v>, A) one has to replace the minus sign in the cosine function by a plus sign. Further- 
more one has to treat the terms where u + \ = N/2 — 1 separately because one state lies at 
the Fermi level and has a different cell width. 



B.2.2 The net integrated density 

n-1 

If we occupy all Wilson states below the Fermi level then we obtain FJ cJ,$o with n = N/2 

u=0 

and $o the vacuum state. This state is not really the free electron ground state. To obtain 
the latter we have to occupy the states (x) and ip^ (x) . They don't interact with the 
impurity but they are occupied. Therefore we define as quasi-vacuum $q the state in which 

all states (fk (x) and (p k (x) with k < /cp = 1 are occupied. Then the ground state is = 

n-1 

\~[ cj,$g. This state has a constant electron density in real space. 

v=0 

n-1 

In the presence of the impurity the new ground state ^ n ew = II a l^'o mus t also contain 

i=0 

this quasi vacuum, i.e., the non- interacting states must be occupied up to the Fermi level. 
Since the inert states are occupied in and ^ new they cancel out when one calculates 
the change in the electron density. The net density of the new state ^ new is the difference 
between p (^ new ) and p (^o)- Since the inert states cancel out one can ignore their existence 
during this calculation. 
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